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ABSTRACT 

Tracer-diffusion of small molecules through dense systems of chain polymers is studied 
within an athermal lattice model, where hard core interactions are taken into account by 
means of the site exclusion principle. An approximate mapping of this problem onto dynamic 
percolation theory is proposed. This method is shown to yield quantitative results for the 
tracer correlation factor of the molecules as a function of density and chain length provided 
the non-Poisson character of temporal renewals in the disorder configurations is properly 
taken into account. 



1. INTRODUCTION 



Atomic charge or mass transport processes in condensed systems often take place in a 
dynamically disordered host medium, whose microscopic structure fluctuates on a time scale 
of the order of the atomic hopping time. An example of current interest in materials science 
are polymer ion conductors.0 These are solutions of ionic salts in a polar polymer that can 
possess significant ionic conductivities. It is well-known that ionic motions in these materials 
are strongly coupled to motions of polymer chain segments, a situation which may be viewed 
as implying a continuing rearrangement of preferred ionic diffusion pathways through the 
host medium. At the glass transition temperature Tg, large scale segmental motions get 
frozen, suppressing long-range ionic diffusion. Other systems where atoms diffuse in a 
reorganizing host medium include permeation of small molecules through polymer filmsil 
or ionic motions through protein channels passing biological cell membranes.^ 

Important progress in calculating the diffusion coefficient of a random walker in a dy- 
namically changing environment emerged from dynamic percolation theory (DPT) and its 
generalizations. In its original form due to Druger et al.i'i one considers the random walk 
in a bond percolation model, where configurations of open and blocked bonds are ran- 
domly renewed at a given rate A. An important outcome of this model is the fact that the 
frequency-dependent diffusivity D{—iuj, A) can be obtained by analytic continuation of the 
diffusivity Do{—iuj) = D{—iuj,0) in the absence of renewals 



irrespective of the precise form of the function Dq{s) to be derived from a system with only 
static disorder .i The same result (|1]) was independently obtained by Harrison and Zwanzig 
within effective medium theoryS assuming independent random renewals of individual bonds 
rather than global renewals as in Ref.i, and also by Hilfer and Orbach.i Subsequent work 
on polymer ion conductors was focused on an identification of the central parameter of this 
theory, the renewal rate A, from experimentally observed polymer viscosities^ and more 
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recently from dielectric relaxation spectroscopy.^ In parallel, the theory was generalized 
considerably to off-lattice hopping|il, spatially correlated renewal^, cases with distinct 
kinds of migration steps^l and, in particular, to non-Poisson renewal processes characterized 
by some waiting-time distribution ■ip{t) .0 In that case the zero-frequency diffusivity in c? = 3 
dimensions is given by 

_ l /-rft^(t)(r^(t))o 

6 f,-dttm ' 

where {r'^{t))o is the mean-square displacement of the random walker in a frozen 
environment.0 Note that in the case ip{t) = Aexp(— At) equation (^) reduces to the zero- 
frequency limit of (|lD. 

While DPT or, generally, dynamic disorder hopping theory0 was developed as a frame- 
work for diffusion of small guest molecules in a fluctuating disordered host environment, it 
was also recognized that the basic idea underlying these dynamically-disordered hopping 
models can provide an approximation to many-particle effects in transport processes in in- 
teracting lattice gases.Ei A (point-like) tracer particle in an interacting lattice gas can hop 
to a neighbouring site provided the other particles have arranged such that this attempted 
site is vacant and that energetic conditions for the hop are fulfilled. The fact that the 
timescales for the changing environment and the tracer motion are interconnected offers a 
way to establish an effective dynamic bond percolation model for the tracer, involving a time 
constant A. A can be determined either self-consistently or by an ansatz based on the lattice 
coordination number. A many-particle effective medium theory for diffusion in interacting 
lattice gases emerges in this way.El 

Besides these investigations for a lattice gas of point particles it seems that DPT-theories, 
although motivated by processes in polymer electrolytes, have never been tested quantita- 
tively in the context of statistical polymer models. While the renewal processes associated 
with a system of point particles are sufficiently characterized by a single rate constant A 
entering equation ([I|),0 we expect this equation to fail for the problem of diffusion through 
a polymer network because of the inherent distribution of relaxation times characterizing 



the chain motion. This should result at least in a more complicated form of the waiting- 
time distribution tp{t). The question now is how far equation can describe the diffusion 
coefficient, when tplt) is defined in a suitable way in terms of the actual dynamics of the 
polymer network. 

To elucidate this question, we investigate in this paper an athermal lattice model, defined 
in Section 2, which consists of lattice chains with varying density and chain length and a 
sufficiently dilute system of point particles. Both chains and point particles undergo diffusion 
via elementary stochastic moves. This model is a special case of a more general lattice 
model of chains and point particles with specific interactions, used previously to describe 
the influence of temperature, pressure and salt-content on diffusion and network relaxation 
properties of polymer electrolytes. 

In the present work, we first obtain diffusion coefficients from dynamic Monte Carlo 
simulation of our model. These results serve as a reference with respect to the subsequent 
approximation method based on dynamic percolation theory. To implement this theory, we 
determine by simulation i) the waiting-time distribution ipit), which we define in terms of the 
occupational correlation function of a site next to a fixed point particle and ii) the mean- 
square displacement (r^(t))o of point particles for static disorder (frozen chains). These 
steps are computationally much less demanding than the full simulation. Comparison of 
both methods via equation @ provides a sensitive test for the applicability of DP-theories 
to diffusion in a fluctuating polymer host. We find excellent agreement between the tracer 
correlation-factors as a function of density and chain length, as obtained from those two 
methods. Temporal correlations, reflected in the non-exponential character of ip{t), are 
found to be crucial in this analysis.^ 

In Section 3 we specialise to a chain length r = 1 which corresponds to a system of 
point-particles only, before we present in Section 4 our full analysis for chains up to a length 
r = 20. Some further conclusions are drawn in Section 5. 
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2. SIMULATION METHOD AND IMPLEMENTATION OF THE DYNAMIC 

PERCOLATION CONCEPT 



Consider a system of lattice chains on a three-dimensional simple cubic lattice of spacing 
a. The chains are made of beads, assigned to lattice sites, and linearly connected via 
nearest-neighbour bonds. Apart from site exclusion, which mimics a hard-core repulsion, 
no explicit interactions between beads are assumed. For chains each with r beads in 
a box of linear size La, the concentration of occupied lattice sites is simply given by c = 
Nr/L^. In addition, our system contains point-like tracer particles, again subjected to site 
exclusion, with a concentration q <^ 1 sufficiently small so that correlations among them are 
negligible. Most of our simulations were carried out with L = 10, r = 1 to 20, q = 10~^, and 
periodic boundary conditions are employed. After preparation of the system with the desired 
number of chains, equilibration and the subsequent dynamics at equilibrium are based on the 
generalized Verdier-Stockmayer algorithm, which employs end-bond motions, kink-jumps 
and crankshaft rotationsEl^. Point-particles individually perform nearest-neighbour hops. 
In the special case r = 2 (moving dimers) only the end-bond motion is active, which then 
is a 90-degree rotation of the dimer about one of its end-points. As usual, introducing 



we can obtain the diffusion coefficient of point-particles, D = lims_>o+ -D(s), from their 
simulated mean-square displacement (r^(t)). To separate the average effect of blocking, 
contained in a factor 1 — c, one introduces the tracer correlation factor /(c) < 1 according 
to 



where D^^'^ = Ta^ denotes the diffusion coefficient for infinite dilution (c — > 0), with T the 
bare hopping rate. 

Our aim is now to map the complete system dynamics onto a disordered single-particle 
model, where disorder configurations are globally renewed according to some appropriate 
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waiting-time distribution ipit). In order to test the validity of this idea against full simula- 
tions, we have to extract the input quantities to equation (Q), (r^(t))o and ip(t), from our 
polymer model. While (r^(t))o can be obtained in a straightforward manner from separate 
simulations with frozen chains, determination of tp{t) requires more explanation. As indi- 
cated already in the Introduction, we propose to determine ipit) from the local occupational 
correlation function (nj(t)nj(O)), where nj(t) is the occupation by chain beads of a site i ad- 
jacent to a fixed tracer position. To simulate {ni{t)ni{0)) chains were first equilibrated while 
keeping the tracer fixed. Such a procedure is perfectly in the spirit of dynamic percolation 
theory based on renewals as "seen" by the tracer in its immediate neighbourhood. Let us 
introduce the probability with t > that there is no renewal within the interval [0,t] 
when the foregoing renewal took place at an arbitrary time to < 0. Following Ref.0, 

= 1- [ dt'(f){t') (5) 
Jo 

withi 

where = (0(0))^^ = dttip(t) denotes the mean renewal time. Now, we argue that 
with probability ^{t) the occupation at site i (next to a fixed tracer) does not change within 
[0,t\ so that the stochastic variable ni{t) (with possible values or 1) preserves its initial 
value, riiit) = nj(0) and nj(t)nj(0) = (?2j(0))^ = ^^(O). Conversely, with probability 1 — $(t), 
one or more renewals occur within [0, if:]. Then, since configurations are randomly reassigned, 
ni{t) can be replaced by its average, c. Hence, in this case, nj(t)nj(0) = cnj(O). Averaging 
in addition over the initial occupation ^^(0), we obtain for the correlation function 

(n,(tK(0)) = c[$(t) + c(l-$(t))]. (7) 

This can be rewritten as 

(n,(t)n,(0)) -c^ 
- ^(1^) ' 



consistent with the requirements $(0) = 1 and — as t ^ oo. Combination of (|5|), 
and (H) yields 

= [\c{l-c)Y'^^{n,{t)nm) (9) 

After insertion into (§) the prefactors drop out. Equation completes the implementation 
of DP-theory to our many-particle model. In the next sections we test the performance of 
this approximation scheme to a simple hard-core lattice gas and to a polymeric system. 
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3. HARD CORE LATTICE GAS 



As a first application let us briefly examine the special case of non-connected beads, 
r = 1, which is identical to the conventional hard-core lattice gas of point-particles. The 
tracer correlation factor /(c) in that case is known to a high degree of accuracy via dynamic 
pair approximations&ii, giving 

rf X ^ 1 + (cos 9) 

' 1- [(3-2c)/(2-c)](cose)' ^ ^ 

and through simulations.il In Ref. which becomes exact as c ^ 1, the quantity (cos0) 
characterizes the average directional change in two consecutive steps of the tracer due to the 
presence of one vacancy. For a simple cubic lattice, (cos0) ~ —0.209. An effective-medium 
approximation to /(c) was obtained recentlyEl from dynamic percolation theory using the 
Harrison Zwanzig approach.0 

In what follows we apply the approach outlined in Section 2 to the same problem. {r^{t))o 
is obtained from simulating a single mobile particle in the frozen conflguration of the back- 
ground particles, ^(t) is deduced from Eq. (^ where {ni(t)ni{0)) is obtained from a short 
time simulation of a lattice gas with one flxed tracer particle, as described above. These 
simulations were carried out within a cubic box of length L = 10 and periodic boundary 
conditions. Note that collective properties of a hard-core lattice gas with symmetric tran- 
sition rates show a relaxational behaviour independent of concentration0ii. The function 
$(t) as determined from (|^) is therefore c-independent and thus can be determined from 
single particle random walk theory. Within that framework $(t) can be interpreted as time- 
dependent probability of return of a single random walker to site i, taking into account that 
one site adjacent to i is blocked by a flxed tracer. In Appendix A we briefly indicate how 
^(t) can be calculated exactly or how one can generate efficient analytic approximations. 

Results for ^(t) obtained both from Monte Carlo simulation and from these approxima- 
tions are plotted in Fig. 1. As seen from the figure, the main decay of $(t) at short times 
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is fairly well represented by an exponential with decay rate 

Ao = -(rf$/rft)i=o = (5/6)r (11) 

The factor 5/6 simply arises from the fact that one of the six bonds in the simple cubic lattice 
connected to site i is blocked by the tracer. The actual decay of ^{t) is approached gradually 
by continued fraction approximants of increasing order A^, which were derived according to 
Appendix A. The asymptotic decay of at long times is governed by diffusion, giving 
^(t) oc (rt)"'^/^, which, however, cannot be accounted for by a finite continued fraction. 
For the purpose of practically evaluating (0) we find that it is sufficient to approximate 
in terms of a superposition of three exponentials. Fig. 2 shows the c-dependent tracer 
correlation factor obtained in this way. The agreement of data points from our DPT with the 
full curve representing the dynamic pair approximation Eq. ( p!OD is quite satisfactory. For 
completeness we also included Monte Carlo data for the full hard core lattice gas. The DPT 
result for c = 1 with value / ^ 0.6802 was obtained analytically, see Appendix A, whereas the 
exact value is /(I) = (1 + (cosB))/(l — (cos0)) ~ 0.654. Also shown are diffusion constants 
calculated from the effective medium approximation as described in Ref.0, which is based 
on only one time constant for renewal events. With A = Aq as given by (|n|) this theory 
yields the dashed curve which deviates notably from Eq. ( p!ol) in the high-concentration 
regime. (On the other hand, merely fitting A to the exact value of /(I), giving A ^ 0.62 P, 
turns out to give very good agreement with Eq. ([10|) in the whole concentration range.) 

These results confirm the conclusion in Ref.lll concerning the applicability of dynamic 
percolation theory to many-particle systems and at the same time indicate that the theory 
significantly improves when the non-Poisson character of renewal processes is taken into 
account. For the problems in the next section this last aspect will become much more 
important. 
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4. LATTICE POLYMERS 



For the hard-core lattice gas (r = 1) the procedure described obviously bears no com- 
putational advantage over existing methods. The situation changes, however, when we go 
over to r > 1. The correlation factor /(c) now depends on r and no analytic approximation 
equivalent to Eq. (|r^) is available for this case. At the same time full simulations of the 
diffusional dynamics become more demanding because of the internal degrees of freedom of 
the host molecules and the larger statistical errors connected with the small concentration 
of tracer particles. Here our approximate DBT-based computational scheme is potentially 
useful. In this section we examine the performance of this scheme. Fig. 3 summarizes our 
MC simulation results, again represented in terms of the correlation factor /(c), Eq. (^, for 
different chain lengths r up to r = 20. The full lines are fits to the simple functional form 
/(c;r) = (1 — a{r)c)/{l — /3(r)c) with fit parameter a, (3 that depend on the chain lengths 
r. These results will be used as a basis for assessing the performance of the approximate 
DP-based approach. As discussed above, this approach is based on evaluating the waiting 
time distribution ip(t) according to Eq. (|]) and the mean-square displacement (r^(t))o of 
a tracer in the presence of a frozen solvent. Figure 4 shows typical results for the function 
$(t), see Eq. (|]) obtained for chains with length r = 10 for several concentrations. While 
our simulation results for (r^(t))o are shown in Fig. 5 for the same r- values as in Fig. 4. 
Substitution of these results for into Eqs. (2) and (4) yields our DP-approximation for the 
correlation factor /(c), which is shown in Fig. 6 together with the 'exact' MC results of 
Fig. 3. (For clarity, only the full lines from Fig. 3, representing the fitted data as discussed 
above, appear in Fig. 6). Evidently, the DP-approximation agrees very well with the full 
simulation for all r. Some further observations are noteworthy: 

(a) From the MC simulation results for /(c;r) (Fig. 3) an interesting picture emerges 
concerning the effect of host connectivity on the tracer diffusion. The most prominent effect 
is the special behavior of chains with r = 2 (see below). Focusing first on chains with r > 5 
we see that for c^0.8@ /(c) is larger than in the hard core lattice gas (Sec. 3), showing that 
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in this density range chain connectivity facihtates diffusion of tracer particles. For larger 
concentrations, however, /(c) drops markedly. 

(b) These findings for r^5 contrast to the behavior found for r = 2 which is a special case 
concerning the allowed elementary moves (see Sect. 2). For all concentrations c considered, 
/(c) now remains larger than 0.9, see Fig. 3. Dimers therefore induce only minor backward 
correlations in the tracer motion. Intuitively, from the point of view of the tracer, only one 
monomer of the dimer molecule effectively suppresses tracer forward motion by a nearest- 
neighbor hop, while the second monomer is shielded. 

(c) The function $(t) is seen in Fig. 4 to decay in a highly non-exponential fashion, 
indicating the importance of temporal correlations in the associated renewal processes. Fur- 
thermore, for small concentrations the relaxation first becomes faster as c increases, which 
reflects enhanced fluctuations after the onset of overlap between chains, whereas in the high- 
density system c = 0.8 the decay is markedly slowed down. For dimers (r = 2) we have 
found that ^{t) decays even somewhat faster than in the case r — 1 and is only weakly 
c-dependent. 

(d) The mean square tracer displacement (r^(i))o in a frozen host, plotted in Fig. 5, 
shows a cross-over from diffusive behavior (r^(i))o ~ to a locahzed random walk (r^(i))o — > 
constant as t — > oo as expected for a percolative network. It is expected that this crossover 
takes place at some critical concentration Ccrit- A precise determination of the percolation 
threshold Ccrit('") for walks through a frozen network of chains of length r is beyond the scope 
of this article, yet rough estimates are presented in Appendix B for 2 and 3-dimensional 
systems. In d = 3 dimensions, Ccrit appears to increase with r, indicating again that for 
given c the frozen chains are less prohibitive to tracer diffusion than a frozen background of 
independent monomers. For example, for r = 10, the concentration c = 0.8 clearly exceeds 
Ccrit('" — 10), (see Fig. 5 and the estimates in Appendix B). 

(e) As already noted the DP-approximation agrees very well with the full simulation for 
all r. This remains true even in the special case r = 2 and has the advantage of saving up 
to about one order of magnitude in computing time. 
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5. SUMMARY AND CONCLUSION 



A method has been proposed how to map particle diffusion through a fluctuating network 
of polymer chains onto dynamic percolation theory (DPT). As input quantities this theory 
requires the particles' mean-square displacement (r^(t))o in the frozen network and the 
waiting time distribution ip{t) for network renewals. We proposed to relate ip{t) to the 
occupational correlation function of a site next to the fixed tracer particle so that it reflects 
the temporal distribution of pathway openings seen by the fixed tracer. In contrast to 
the standard hard core lattice gas, ip{t) decays in a highly non-exponential fashion when 
longer chains are considered. This feature of the fluctuating network appears to be crucial 
in implementing dynamic percolation theory to chain systems. When properly taken into 
account, the DP-model gives quite accurate results for the tracer correlation factor in its 
variation with concentration c and chain length r. We have verified this by comparing the 
results of DPT with Monte Carlo simulations of the complete system dynamics. 

Our studies so far are limited to an athermal system. The theory in that case was found 
to save about one order of magnitude in computing time relative to full simulations. Under 
this aspect it would be very interesting to extend these studies by applying DP-theory 
with non-Poisson renewals to interacting systems and associated questions of dispersive 
transport. 
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APPENDIX A: HARD CORE LATTICE GAS: SUPPLEMENTARY RESULTS 

Regarding the hard-core lattice gas, the quantity $(t) given by (H) is equal to the 
probability of return to site i after time t in the case of a single random walker. Here site 
i is next to a site blocked by the tracer. In calculating $(t), the effect of the blocked site 
can be taken into account by standard defect matrix theory.0 This yields an expression for 
in terms of the solution of a 3 x 3 matrix equation, whose coefficients are determined 
by unperturbed lattice Green functions up to third-neighbor distances. The calculation is 
straightforward and will not be reproduced here. 

In the present context it is sufficient to obtain an accurate approximation for only 
until it decays to about 10^^. Computationally it is then advantageous to represent its 
Laplace-transform <l>(s) as a continued fraction of the type <l>(s) = ao(s+6i— ai(s+. . . )~^)~^, 
generated by a short time expansion of Time derivatives = 0) with < n < 

2N — 1 are easily obtained by enumerating closed paths of the walker which avoid the blocked 
site. Specifically, we use 

m = f]^^-^e-^' (12) 

n=0 

where is the probability of return to the origin after n steps. 

At stage the continued fraction is terminated such that $(0) agrees with the exact 
result from defect matrix theory (see above) which, at s = 0, is determined in terms of 
Watson-type integrals. Fig. 1 contains a plot of the A^-th order approximants for $(t) up 
to = 6. For = 6 the simulations are accurately represented up to t ~ 15. 

Finally we comment on the limit c — > 1, where the correlation factor from DPT can 
be evaluated analytically. Obviously, the mean-square displacement (r^(t))o of a tracer in 
a frozen lattice is determined in that limit by successive exchanges with one neighboring 
vacancy. This gives (r^(t))o = (1 — exp(— 2rt))/2, where F is the jump frequency. From 
(§) together with the above-mentioned results for $(t) we obtain /(I) = 0.6802 as given in 
Section 3. 
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APPENDIX B: PERCOLATION IN A FROZEN NETWORK 



Calculation of the diffusion constant D from requires knowledge of the mean square 
displacement (r^(t))o of a tracer particle in a frozen network on time scales of the order of 
the decay time of il){t). In this Appendix we estimate the critical concentration Ccriti'f') for 
percolation of a monomer particle through a frozen network of chains, which distinguishes 
diffusive from localised behavior of (r^(t))o at long times. To our knowledge, this problem 
of correlated percolation has not been investigated before for general chainlength r. We do 
not, however, attempt any precise determination of Ccrit{f^)', rather we like to point out some 
major qualitative trends in Ccrit{f^) as a function of r in c? = 3 and d = 2 dimensions. 

In d = 3, equilibrated chain configurations were prepared by the same algorithm as 
described in Section 2, while in ci = 2 we used the algorithm by Siepmann et al.il For 
systems of varying size L, increased in steps AL, we determined the probabilities P(r, c, L) 
of occurrence of a spanning cluster of vacant nearest-neighbor sites. For fixed r, we obtained 
points of intersection of successive curves P(r, c, L) and P(r, c,L~ AL) versus c, which give 
successive approximations for Ccrit{r)^ Estimates for Ccrit{r) were deduced from calculations 
with L-values up to a maximum Lmax, to be adapted to r. For example, in c? = 2 we went up 
to Lmax = 200 for r = 20, whereas L^ax = 80 for frozen dimers r = 2. These values appeared 
sufficient to achieve reasonable convergence. In the limit r = 1 we find Ccru — 0.41 for d = 2 
and Ccrit = 0.69 for d = 3, which reasonably agree with values 1 — pc given by the well- 
known thresholds pc for site percolation on square and simple cubic lattices, respectively. 
The qualitative r-dependence of Ccrit{r) is shown in Fig. 7. In d = 3, Ccritij") monotonously 
increases with r which we interpret as a reduction of blocking of open pathways through 
the connectivity of chains. The most pronounced increase occurs already when going from 
monomers to r = 2. \n d = 2 this argument again applies to the step from r = 1 to 
r = 2, but blocking becomes more effective for longer chains so that Ccrui'f') decreases, with 
Ccrit{r) < Ccriti^) foT r^l2. For long chains, especially in d = 2, it might be interesting 
to study (r^(t))o near criticality on different length scales above and below the radius of 
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gyration of chains, but this goes beyond the scope of this work. 

Generally, for sufficiently long times and small |c — Ccritl such that the correlation length 
becomes larger than the size of chains, one expects (r^(t))o to follow the standard scaling 
forms for non-correlated percolation, which imply associated scaling forms for the frequency- 
dependent diffusivity Dq^—iuj) as — ^> 0. In applications of DP-theory where the analytic 
continuation rule (|lD holds at least for small A, one can immediately predict the asymptotic 
forms of the long-time diffusion constant D = -Do (A) of a walker in the presence of slow 
(A — > 0), Poisson-type network renewals.0 The resulting scaling expressions for D are 
straightforward to write down from the corresponding expressions for Do{—iuj) given in 
Ref.@. Some notable special cases for nonzero A are 

{A , C Ccrit, 

(13) 
A I C ^crit I , C ^ ^crit 

where ly and f3 are the conventional static percolation exponents for the correlation length 
and the order parameter, respectively, and k is the dynamic critical exponent for anomalous 
diffusion in a percolation system at criticality. 
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FIGURE CAPTIONS 



Fig. 1 Semilogarithmic plot of the function (see equation 8) for the hard core lattice 
gas. Monte Carlo data are obtained by simulating the probability of return of a single 
random walker starting next to a blocked site. 

Dashed-dotted line: single exponential approximation as determined by the initial 
slope (see equation 11). The other thin lines represent continued fraction approximants 
up to order = 6 (see Appendix A). 

Fig. 2 Tracer correlation factor /(c) of the hard core lattice gas against concentration, ob- 
tained by different methods (see text). 

Fig. 3 Simulated tracer correlation factor /(c) for different chain lengths r = 2,5, 10 and 20. 
Full lines refer to the fit function / ^ (1 — ac)/(l — /3c), where for r = 5, 10 and 
20 the fit parameters are /? = 1 and a = 1.057, 1.062 and 1.071, respectively. This 
implies that the diffusion constant (^) is approximately linear in c, -D(c) ~ -Do(l — «c). 
One the other hand, for r = 2 we find a = 0.391 and f3 = 0.318. For comparison we 
also show simulation data for r = 1 (hard core lattice gas) together with Eq. (p!OD 
(dashed-dotted line). 

Fig. 4 $(t) for chains of length r = 10 for three different concentrations c = 0.8, 0.1 and 0.4 
(from above). Also shown is the short time behavior of $(it) for r = 1 (dashed-dotted 
line), reproduced from Fig. 1. 

Fig. 5 Mean-square displacement {r'^{t))o of walkers in a frozen chain network. Parameters 
are as in Fig. 4. For c = 0.4 a comparison is made with the case r = 1. 

Fig. 6 Comparison of tracer correlation factors from DP-theory for chains of different lengths 
(data points) with results from full simulations. Full lines represent fit functions for 
the simulation data, reproduced from Fig. 3. 
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Fig. 7 Estimates for the critical concentration Ccrit{r) versus r for percolation of a monomer 
particle through a frozen network of chains of length r on simple cubic {d — 3) and a 
square {d = 2) lattices. 
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